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Abstract 

We discuss nonperturbative contributions to the 3-dimensional one-loop ef- 
fective potential of the electroweak theory at high temperatures in the framework 
of the stochastic vacuum model. It assumes a gauge-field background with Gaus- 
sian correlations which leads to confinement. The instability of < F 2 >= in 
Yang-Mills theory appears for small Higgs expectation value < 4> 2 > in an IR 
regularized form. The gauge boson propagator obtains a positive momentum- 
dependent "diamagnetic" effective (mass) 2 due to confinement effects and a neg- 
ative one due to "paramagnetic" spin-spin interactions which are related to the 
< F 2 >= instability. Numerical evaluation of an approximate effective poten- 
tial containing these masses shows qualitatively the fading away of the first-order 
phase transition with increasing Higgs mass which was observed in lattice cal- 
culations. The crossover point can be roughly determined postulating that the 
effective 4> A and d> 2 terms vanish there. 
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1 Introduction 



The phase diagram of the electroweak standard model at high temperatures T of about 
100 GeV can be well described by the 3-dimensional effective Higgs field gauge theory 
which is obtained perturbatively by "integrating out" heavy degrees of freedom, i.e. 
by matching a set of amplitudes in 4 and 3 dimensions [p], [|: 



C 3 = \f%F« v + {D^) + {D^) + + A 3 (</> + 0) 2 . (1.1) 

Here higher derivative terms are neglected. This is adequate to order g^, where g w is the 
weak coupling constant in d = 4. The SU(2) Yang-Mills action ^F 2 and the covariant 
derivative contain a gauge coupling g\ proportional to g%T. The U(l) part of the 
standard theory is neglected here since it is not important for our considerations, but it 
could be easily included. Furthermore, and A3 are the effective T-dependent mass 
and coupling constant of the dimensionally reduced theory, respectively. 

The dimensionless ratio g1/m\R with some typical physical infrared (IR) scale miR, 
is not small in general. Indeed, the gauge coupling in a Wilson-type effective action 
increases strongly with the IR cut-off going to zero ||. Thus the further use of the 
perturbation theory based upon the Lagrangian (|1.1|) is very dangerous (like in QCD) 
and it is only consistent if the scale < + > is not too small. Besides this, one also has 
to check if the quantity which one tries to calculate perturbatively is not dominated 
by nonperturbative contributions, in particular if the former turns out to be small, as 
in the case of the effective potential to be discussed. 



Since the fermions are integrated out already, the action (1.1) can be easily put on 
a lattice [|, ||. The critical temperature T c of the phase transition, v(T c )/T c with v(T c ) 
the Higgs field minimum, the latent heat, and the surface tension can be calculated 
("measured") with good accuracy. The dimensionless quantities 

V=4> (1-2) 

03 03 

span the phase diagram; y is connected to (T — T c ) and x describes the strength of 
the phase transition and depends on the Higgs mass. For x < 0.03 — 0.04 the first 
order phase transition obtained in two-loop perturbation theory compares very well 
with lattice calculations 0] of the above-mentioned quantities. 

For larger x values there are deviations which become stronger for quantities of de- 
creasing perturbative order in g w . In the above list, the surface tension shows the worst 
failure of a perturbative calculation ||. For x < 0.11 (corresponding to m# ~ raw) 
lattice calculations show || that the phase transition fades away, there is a crossover. 
On the other hand, a two-loop perturbative treatment |], 0] gives a first-order phase 
transition. Perturbation theory can be formally applied since g%/g w v(T) is still rather 
small in this x-range; but the effective potential becomes very shallow and v(T c )/T c is 
quite small, so that the phase transition is very weakly first order. We conclude that in 
such a situation nonperturbative contributions to the effective potential are essential. 

In this paper we develop a physical picture for such nonperturbative contribu- 
tions. We will present analytic indications for an instability of the vacuum of the 
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3-dimensional theory (|1.1|) at < F >— in case of small Higgs vevs. We then will 
develop a up to now qualitative model for the < F 2 >^ vacuum contribution to the 
electroweak potential. This will lead to an understanding of the crossover. 

Lattice studies reveal that in the high temperature phase of the effective theory 
( [LID , i.e. at values of m| corresponding to T > T c , one has confinement phenomena 
like in QCD. There is a linearly rising part in the potential between static sources || 
and there exist PF-ball states almost identical to those in pure SU(2) YM theory |§, 
which fits nicely together with the observation that the string tension is approximately 
the same in both theories. Also three-dimensional correlation masses of + ,l~,2 + 
bound states of Higgses have been "measured" ||. They can be described || in a 



simple relativistic bound state model analogous to the one of Simonov |T(| in QCD 
whose main ingredient is an area law for the Wegner- Wilson loop. 

This area law follows naturally from the stochastic vacuum model of Dosch and 
Simonov JTT], [12| which is quite successful in QCD. It postulates the following gauge 



and Lorentz covariant correlation function between field strengths = F^ u T a at 
different positions x and x' 



X 



(S^ua - 5 iut 5 v x)kD{z I a ) 



+2-^ (d u (zj, x - z x 5, a ) - (/i <- u^D^/a 2 )] (1.3) 

to dominate the cumulant expansion (see eq. (|3.4f) ). Here < g 2 F 2 > is the usual 
x-independent condensate, and D and D\ are form factors, containing a correlation 
length a normalized such that D(0) = -Di(O) = 1. The way it is used in this paper 
( P-.3|) only contains the nonperturbative correlation in the nontrivial QCD type vacuum 
(we skip the index NP in the following). Thus in the background field formalism 
with a decomposition A® = ^"(background) _|_ a a(quantum) ^ describe a correlated gauge 

field background. To obtain the nonperturbative correlation in the continuum limit 
from lattice calculations the (diverging) perturbative part has to be subtracted (or 
suppressed by "cooling"). The quantum field a® propagates and couples according 
to the usual perturbation theory. We expect that the propagation of the soft quanta 
should be suppressed by the background and this will indeed turn out to be the case. 
Furthermore 



F* v {x,x Q ) = Pexp(ig / dx x A x (x)) 

l Jx 



ab 

FlM (1.4) 



is the field strength tensor parallel-transported to a fixed reference point Xq by a 
Schwinger string in the adjoint representation. The integration path is fixed to be 
a straight line. But even then the l.h.s. of ( |1.3| ) generically depends on xq whereas 
the r.h.s. only depends on z = x' — x. Neglect of the xq dependence is an assumption 



12| , which is fulfilled in an appropriate range of xq where the cumulant expansion con- 
verges rapidly: The Wegner- Wilson loop vev is Xq independent and if it is evaluated 
with stochastic correlations also these have to be x independent. With a linear path 
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in ( |1.4| ) and with xq on the line x — x' the ansatz ( |1.3| ) has been tested in 4-dimensional 
QCD lattice calculations |13| (and is being also probed in recent 3-dimensional lattice 
gauge theory investigations fl4j0. 

In a coordinate (Fock-Schwinger) gauge with reference point at x$ the parallel 
transport operator in ( |1.4|) , along a straight line, equals 5 ab . Using this gauge in the 



correlator ( |1.3| ) (or more general the nonabelian Stokes theorem) in order to perform a 



cumulant expansion of the vacuum expectation value of the Wegner- Wilson loop with 



the higher cumulants neglected, one easily obtains [11, 12] the area law. The string 
tension is related to the local condensate < g 2 F 2 > and to the correlation length a 
which enters the form factors D and Di (see eq. ( |3.6D ). 

Applying this picture to the three-dimensional Higgs-gauge theory a stochas- 

tic gauge field background would not only be present in the "hot" phase, but more 
generally at small Higgs background fields <fi. It leads to a modified effective poten- 
tial containing nonperturbative effects of the gauge field background. This has been 
already proposed [[OJ some time ago for constant < F 2 >. In the present paper the 
importance of nontrivial form factors describing a correlated gauge field background is 
demonstrated. 

We propose that the nonperturbative dynamics of the theory (1.1) is dominated 
by a fluctuating IR-gauge field background like in pure YM theory but with the Higgs 
field background as a further parameter. Our most important starting point will be 
the instability of a vacuum with < F 2 >= for small background Higgs fields. Such 
an instability we obtain due to paramagnetic interactions related to the spin of the 
'W-gluons". It is similar to the Savvidy instability of QCD. But it is now stated for 
fluctuating fields, and there are no IR singularities. 

Chapter 2 contains a discussion of the order F 2 -term in an effective potential V 
depending on the gauge and Higgs fields and of the instability of a vacuum with 
< F 2 >= 0. In chapter 3 we calculate the 1-loop effective potential in a certain 
approximation. We introduce a momentum dependent positive (mass) 2 m 2 on j of the 
gauge boson due to confinement and a negative momentum dependent (mass) 2 — E re- 
lated to paramagnetic spin-spin interactions. The evaluation of the potential in chapter 
4 allows us to qualitatively discuss the crossover behavior of the hot electroweak theory. 
Appendix A contains a thorough discussion of the gluon propagator in the stochastic 
vacuum background. While the basic idea of how stochastic background fields can lead 
to the formation of condensates is very easy to understand, its implementation in a 
gauge theory is rather involved. Therefore, in Appendix B, we discuss this mechanism 
for the much more transparent case of a simple scalar theory. The reader might wish 
to turn to this Appendix before embarking on the detailed presentation in section 2 
and 3. 



2 Instability of F 2 = for small Higgs field (p 

As a first step let us calculate the term 

< T FF >= J d 3 xV FF (if) (2.1) 



3 



in the 3-dimensional effective action which is linear in the condensate 
< g 2 F 2 >=< giF« v F« v >: 

V FF (<p) =< glF 2 > P{ V 2 ) (2.2) 



where <p = y/2 < >• The potential Vff arises from the stochastic average of that 
term in the ordinary effective action T[A, tp] which contains two powers of F^ u and an 
arbitrary number of covariant derivatives acting on them: 

< T[A, tp] > = < -~g 2 J d d x tr c F M ,(x)n(- J D 2 )F^(x) + ... > (2.3) 

= <r FF >+... 

Here II is the y?-dependent polarization function (divided by g 2 ) and tr c denotes a 
color trace in the adjoint representation. In order to be slightly more general, we shall 
consider a SU (N c ) gauge theory in d dimensions and set N c = 2, d = 3, g = g% only at 
the very end. To proceed, we evaluate the gauge-invariant action ( |2.3j ) for gauge fields 
which satisfy the Fock-Schwinger gauge condition 

(x-x Q )^(x) = (2.4) 

for some arbitrary point Xq. Gauge fields satisfying Q2.4|) can be expressed in terms of 
the corresponding field strength according to 

A a (x) = f dr]r]{x - x ) v F"(x + r}{x - x )) (2.5) 
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As a consequence, there is a one-to-one correspondence between powers of A^ and 
powers of F^, and being interested in F 2 -terms only we may ignore the A^-terms in 
the covariant derivatives of (12. 3D 



r[A,<p] = \g 2 N c J d d xF; v {x)U{-d 2 )F; v {x) + ... (2.6) 



A priori the reference point x in the stochastic correlator ( |1.3|) is unrelated to the 
point xq above which characterizes a specific Fock-Schwinger gauge. In the following 
we shall identify these two points. This has the consequence that the parallel-transport 
operator in eq. (|1.4| ) becomes equal to the unit matrix. Only in this case we are dealing 
with the correlator of two local field strength operators. In particular, 



« g 2 F^(x')F* u (x) »=< g F > D eS [ - ] (2.7) 



with the abbreviation 

°- (3 - " c (3 + <> - ^ (fl + - ^ (?) < 28 > 

Now it is straightforward to compute <C Tpp 3> by applying (|2~7l ) to (|2~^ ): 

Vff(<p) = \Nc < 9 2 F 2 > J ^- d U(p 2 )D cS (p 2 ) (2.9) 



Here D c g(p 2 ) is the Fourier transform of D e ^{z 2 /a 2 ). Denoting the Fourier transforms 
of D(z 2 /a 2 ) and D\(z 2 /a 2 ) as D(p 2 ) and Di(p 2 ), respectively, it reads 



D cS (p 2 )=KD(p 2 )--(l-K)p 2 D[(p 2 ) 



(2.10) 



Later on in our numerical computations we shall consider the case k = 1 onlyQ and 
use the two form factors 



D« 4 



o-\Al a 



-z 2 /a 2 



Their Fourier-transforms for d = 3 are, respectively, 



D (1 \p 2 ) 

b^ 2 \p 2 ) 



d 3 ze ipz e~ lzl/a 



7r 3/2 a 3 exp 



8vra 3 



1 + a 2 p 



-p 



(2.12) 

(2.13) 
(2.14) 



It remains to compute the polarization function EL We restrict ourselves to a one- 
loop calculation here. The dominant contributions to II come from gauge boson and 
ghost loops. (Higgs loops play a minor role and are considered later on.) Using the 
proper time method, their contribution to the effective action is given by 



1 



-Trln[K] - Trln[-D 2 + 



m 



°° dT f 1 



T 12 



-Tr[e- T ^-Tr[e- T (- D2+m2 >i 



Here 



K 



D 2 5n V + 2igF lw + m 2 5 i 



ub 



(2.15) 



(2.16) 



is the kinetic operator of the gauge bosons in the Feynman-'t Hooft gauge and — D 2 +m 2 
the corresponding one for the ghosts. We identify 



2 1 2 2 



2<6+6> 



(2.17) 



Using the method of Barvinsky and Vilkovisky |TJ| or the corresponding world line 
technique [[HJ it is straightforward to identify the F 2 -terms which are contained in the 
functional traces of (|2.15|) : 



Ii{p 2 ) 



( 4vr )rf/2 



r dTe- m2T T- d ' 2+l 
Jo 

f{p 2 T) - 1 



/(pT) + -(d_2) 



p 2 T 



(2.18) 



1 The D\ form factor does not contribute to the area law (ll], |l2) 
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with the familiar "second-order form factor'' 



f(p 2 T) = / da exp[— ail — a)p 2 T] 
Jo 



(2.19) 



Performing the proper time integral yields IT 
(d - 2)r(l - d/2) 1 



n G b 2 ) 



( 4vr )d/2 

4r(2 - d/2) 
(4vr) rf / 2 



JT ^0 



da 



da 



I1 G + Hf with 

" m 2 + a(l - a)p 2 f 2 - 1 

d/2-2 



m 



d-2 



m + a(l — a)p 



^2.20) 



The contribution lip originates from the nonminimal coupling of the gauge boson 
fluctuations to the background (the i^-term in K^), while Uq arises from the D 2 - 
term in K^ v and from the ghosts. Using ( |2.20|) in ( (2.9|) for d = 3 and N c = 2 we obtain 
the desired expression for the potential: 



V FF (<f) 



da 



< glF 2 > 



m 2 + a(l — a)p* 



dpp 2 D eS (p 2 

1/2 1 



(2.21) 



2p 2 



(m + a(l — a)p ) 



2N1/2 



m 



The momentum integration in (p.21 ) converges both for p — > oo and for p — > 0. In fact, 
Vpi? is IR-finite even in the limit m = g-^/2 — >• of massless gauge bosons and ghosts. 
We emphasize that this IR finiteness has nothing to do with the dynamical generation 
of a mass for these particles (see below), but rather with the "inclusive" nature of the 
quantity at hand. For d = 3 and m = the polarization function 



n( P 2 



2 V 16 ; 



(2.22) 



is singular at p — 0, of course, but the volume element d 3 p = 4irp 2 dp in (|2.9| ) has the 
effect of completely suppressing the contribution of p = to Vff{<4>)- (D e &(0) is finite 
and nonvanishing. ) 

In eq. ( p.22|) the "1" comes from Uf while the "—1/16" is due to He- We observe 
that the nonminimal F^ u interaction of the gauge boson fluctuations is the dominant 
effect. The contribution from the minimal couplings of gauge bosons and ghosts is 
smaller by more than one order of magnitude and has the opposite sign. This situation 
persists also for m > 0. 

It is important to note that the one- loop contribution Vff{^p) is negative. It has to 
be added to the stochastic average of the positive tree-level term, 



« ~F%F% »=< 9tF 2 > /Ag 2 . 

If it is large enough so that the sum is still negative, this indicates that the vacuum 
with a vanishing gauge field condensate is unstable and that the true ground state of 
the system is characterized by a nonvanishing value of < F 2 >.Q 

2 See ref. pif for a related discussion in the exact renormalization group framework. 
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At the linearized level, the dynamics of the gauge boson fluctuations a M (x) is gov- 
erned by the Lagrangian 

C = la^ \-D 2 5^ + 2igF tlv + m 2 cU a u + 0{a\) (2.23) 



2 



The coupling a^F^ y a v is the spin-1 analogue of ipa ■ Bip for spin-^ particles. It gives 
rise to a "paramagnetic" interaction of the gauge bosons, while the terms resulting 
from the minimal substitution, a ll D 2 a fl lead to the more conventional "diamagnetic" 
behavior ||18|| . Thus, using the terminology of solid state physics, we can say that the 



instability found above results from a dominance of paramagnetic over diamagnetic 
effects. 

It is interesting to note that the possible existence of such an instability is closely 
related to the fact that the gauge bosons carry spin 1 and have a (tree level) Lande 
factor ql = 2. This is most easily seen if we expose a "colored" particle of spin S — % 
or S — 1 with gL = 2 to a covariantly constant color magnetic field of strength B along 
the z-axis. Its squared one-particle energies are given by |T9[1 

E 2 = m 2 +p 2 z + gB(2n + l) -2gBS z (2.24) 

where n = 0, 1, 2, ... enumerates the Landau levels and S z is the spin projection along 
the z-axis. In eq. ( |2.24|) the terms gB(2n + 1) and —2gBS z stem from the diamag- 



netic and the paramagnetic interactions, respectively. An instability is signalled by an 
imaginary part of E. In the case of fermions we have the two possibilities S z = —1/2 
and S z = +1/2. In the former case, E 2 is strictly positive for all n, but in the latter 
the paramagnetic term precisely cancels the diamagnetic one at n = 0. For m 2 = 
one finds at p z = a marginally stable state at E 2 = 0. For a spin-1 particle one has 
S z = 0, ±1. Hence, for the spin properly aligned to the external magnetic field, the 
diamagnetic energy of the lowest Landau level, +gB, combines with a paramagnetic 
contribution —2gB to yield a negative total magnetic energy — gB. Hence E 2 < for 
m 2 +p1 sufficiently small. (This effect is responsible for the instability discussed in the 
context of the Savvidy vacuum of QCD, for instance.) 

The effective action which we have calculated above describes the interaction of 
gauge boson fluctuations with non-constant external fields of arbitrary momentum. 
Hence their spectrum is certainly not given by ( |2.24j ). However, what continues to be 



true is that for the spin-1 field the paramagnetic interaction overrides the diamagnetic 
one - which would not be possible for Dirac fermions. 

When (p 2 is increased from ip 2 = to larger values, the potential Vff{^) approaches 
zero from below. Thus, starting from < F 2 >^ at ip 2 = 0, one will return to a 
situation with vanishing condensate for ip 2 larger than a certain critical value ip 2 . 

In Appendix B we discuss a scalar toy model and give a simple physical explanation 
of why stochastic background fluctuations tend to trigger the formation of condensates. 
Taking a — > oo in eq. ( [2.21 ) with form factors ( [2.13|) or ([2.14|) leads back to the 



perturbative result ~ g\ < F 2 > /m, implying the 1-loop gauge Z-factor| 

(2.25) 



7 - 1 _ M. 

■^gauge 



7 2 11 

mw Srrigh Smcb 2Am H 



3 There is a misprint in formula (A. 12) of ref. 
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Figure 1: The function G(ma) appearing in Vpp of eq. (2.26) for the two different form 
factors 

with 

m w = m lh = \glv 2 ; ™ 2 Gb = \a 2 ¥ 2 + ^sV 2 + ™\\ ™ 2 h = 3A 3 <^ 2 + m 2 3 

It has an IR singularity at m 2 = \g 2 f 2 = which is removed in our calculation with 
a gauge field correlator. Indeed for a — > oo we obtain the W^-boson plus ghost part of 
V FF as — (23/127r)flf (jF 2 ) as given in eq. ( ggg ). 

The scalar Higgs contributions in the loop can be also taken into account in a 
Barvinsky-Vilkovisky type formula. Higgs fields in the loop are easily included in Vff 
by additional terms in II proportional to (f(—Td 2 ) — l)/(—Td 2 ), and by modified 
prefactors ^fG~ m H T and gfre _m G<> T for the Higgs boson and the Goldstone bosons, 

respectively, instead of (d— 2)e~ m2T . This is in agreement with ( p. 25 ) for a — > oo (iV c = 



2, d — 3). It is obvious that these parts are small compared to the gauge parts. Mixed 
gauge boson-Higgs field graphs do not contribute. 

Evaluation of Vff i n ec L- (|2.9|) for k — 1, i.e. with a pure D-form factor gives 



V FF = - < g\F 2 > aG(ma) (2.26) 

with 



G (ma) = J q °° dpD(m 2 a 2 p 2 ' 



15 _ l\ ( p 
—p — — I I arcsm l n 



,4 pj \ t/p^+A 2 

where D(p 2 ) = a 3 D(m 2 a 2 p 2 ) with p = p/m. For m 2 —>■ this simplifies to 

1 5 r°° - 
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Figure 2: Higgs field interaction with stochastic gauge field background 



The first term inside the curly brackets of eq. ( 2.21|) derives from the F^ v part in K^ v of 



eq. ( [2. 16 ). It contributes 4parcsin(...) in the bracket in (p. 27 ) and dominates Vff- It 



is responsible for the negative sign of Vff which is capable of generating an instability 
at F 2 = 0. 

The numerical evaluation of (|2.27j ) for 



1) D^(k 2 ) = 8vr/(l + k 2 ) 2 and 

2) D^(k 2 ) = n 3 / 2 exp(-k 2 /A) 

corresponding to the form factors (|2.11|) and fl2.12j ), respectively, is shown in fig. |l|. 
The massless limit is given by ( |2.28| ). G(0) of eq. ( |2.27| ) is calculated as 

G (1) (0) ~ 0.1492 (2.29) 

and 

G (2) (0) ~ 0.1322 (2.30) 

respectively. Note that defining a = J °° dzD(z 2 ) one has a = a in case (1) and a = 0.85a 
in case (2); thus aG(0) has a very similar value in both cases if expressed in terms of 
a. 

There is a further contribution to Vff by a Higgs field tree term (fig. |2|). Because of 
the correlator it is similar to a 1-loop term. The gauge-boson coupling to an external 
constant Higgs field is given by {—igd^A a ^){^ + T a ip, Wl — Lp+ u T a &), where <£> qu is the 
quantum fluctuation. The correlator 

« d.A^d.A^x') » (2.31) 

can be conveniently evaluated in the Fock-Schwinger (coordinate) gauge Q2.5|) (like the 
Barvinsky-Vilkovisky form factors jT6[) with xq = (since we are discussing the 
contribution to an effective action with quasilocal terms). Note that the differentiation 
in (|2.31|) should not affect the x, x' in x - Fourier transformation to momentum space 
gives 

VfF = I < £F> > *TT-» JL jT dqq [ (2.32) 
lim (I log <? + "!! + m l) (ZD (*) + ^ < ^ 



p^°\P (p — q) 2 + m 2 J \ \(r] + T]) 2 I (rj + rf) 2 \{rj + rj) 2 
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Figure 3: I q (m 3 a) contributing to Vpp e in (2.33) for the two form factors (2.11) and 
(2.12) indicated by (1) resp. (2) 



The limit lim p ^ (---) equals 8q/(q 2 + m§). After an integration by parts of the second 
term with respect to q, the 77 — 77 integration can be done analytically. The (rj + f\)/1 
and q integrations were done numerically for the special form factors ( 2.1 3| ) and ( [2.14|) . 
This results in 

Vff = ^a 2 I q (m 3 a) < g 2 3 F 2 > (2.33) 

where 

6 ,_ g 4 



I q {m 3 a) = — — / dq 2 N 2 / ??(g) 
(2^)-^ Jo (q z + a z m 3 y 



with some function I v {qa) which can be extracted from Q2.32 ). The function I q (m 3 a 



is plotted in fig. |[ Note that I q {m 3 a) does not exist for m| < 0. I q is a monotonically 
decreasing function of am 3 with I q (0) ~ 0.005 — 0.01. Vpp e is positive and thus supports 
the vanishing of the gauge condensate for increasing ip 2 . 

The potential Vpp of eq. ( |2.26|) and the additional small positive pieces including 



the Vpp e discussed above have to be added to the tree potential < \F 2 >. From eq. 
( |2.26| ) and (|2.33|) we obtain an instability at F 2 = if 



g 2 aG(ma) - m 2 a 2 1 ~ q (m 3 a)/ '18 > 0.25 (2.34) 



An instability at F 2 = requires a correlation mass 1/a < 4G(0)g 3 ~ 0-Qg 3 . 



Lattice studies [J4| indicate 1/a ~ 0.73g 3 ~ w g i ue b a ii/2, as suggested in ref. 



Strictly speaking, our calculation does not give an instability with such a value of 1/a. 
However, keeping in mind the limitations of the present semi-quantitative approach, 
we have to be aware that it might be not very accurate. It is also gauge dependent. In 
perturbative calculations of the electroweak potential one observes that a strong gauge 
dependence of the 1-loop order result indicates that the next loop order contribution is 
needed and that one can arrive in 2-loop order at a quantitative reliable result without 
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much gauge dependence [I], [/J. Thus the gauge dependence of the 1-loop expression is 
a very useful hint. A remaining gauge dependence of the potential [0] together with a 
similar one of the Z-factor in the kinetic term should drop out in the calculation of the 
action of an extremal field configuration such as the critical bubble. But this has not 
really been possible to check up to now even in pure perturbative calculations. Here we 
employed the Feynman-'t Hooft gauge. In our present treatment the use of a specific 
model for the IR sector of the theory could weaken, but certainly not erase the gauge- 
fixing dependence. A 2-loop calculation in the perturbative part is very demanding, 
but would be very interesting. 

A negative F 2 term then leads to a ground state with < F 2 >^ to be stabilized 
by confinement forces in an effective potential V((p 2 , < g 2 F 2 >) as we will argue in the 
following chapters. The gauge boson condensate is determined as the minimum position 
of the "potential" in F 2 . Given the minimum value of this potential, its negative linear 
part at small F 2 allows a rough estimate of the minimum position (e. g. allowing for 
an F A term). Thus the lattice determination of the nonperturbative finite part a 4 g 2 of 



the free energy in ref. [31] in pure Yang-Mills theory (m = 0) may be used to connect 



the F 2 coefficient discussed in section 2 with the gauge condensate and with the string 
tension also known from the lattice. The small negative a 4 reported in [31] then points 
to a negative coefficient of \F 2 with an absolute value much smaller than 1. This 
strengthens our view that we are close to instability with the lattice value of a and 
that only small further (2-loop) contributions are needed to get complete agreement. 
A careful analysis including lattice renormalization is required. 

The instability of < F 2 >= is strengthened with growing g\. If we define the IR 
running of g 2 cS via the prefactor of the F 2 term in the effective action, the negative sign 
of the 1-loop contribution tells us that g 2 aff is indeed larger than g 2 and that it would 
diverge at the border of instability of < F 2 >= 0. In the Wilsonian renormalization- 
group approach of ref. j| g% s (k) at a running IR scale k is the central quantity. Its 
divergence at some k = k c signals confinement like in QCD. It is obtained from a 
differential renormalization-group procedure and certainly goes beyond a 1-loop order 
calculation. One way to look at the effective action (average action) of this approach 
is to supplement the usual background field formalism with a smooth IR cut-off k 



for the momenta of inner propagators [p2j| . The stochastic vacuum with a correlation 



length a (together with the condensate < F 2 >, both in principle to be determined 
dynamically) then constitute some "physical" value for the IR scale k. 



3 The effective action with stochastic 
vacuum correlations 

In the previous section, only those terms in << T[A, if] » were calculated which are 
linear in < g 2 F 2 >. Let us now turn to the much more complicated question of how 
<< T >> depends on < g\F 2 > to all orders. 

The 1-loop effective action with a gauge boson circulating in the loop and with a 
Higgs field plus a correlated gauge field background can be easily written down as a 
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worldline path integral: 



W^<^ » = -Iff / (0) ^[*]exp(- /> 







4 



< tr cL Pexp j- jf dr(-ig) (A M (x)x M ll L + 2F(x))| > 



(3.1) 



The subscripts "c" and "L" refer to the color and Lorentz matrix structure, respectively. 
We use the coordinate gauge based at the center of mass point x |I7| defined by 



x(t) = x + y(T) 



[ T dry(T) =0, Dx = d d x Dy 
Jo 



which amounts to 



^o + y) = J o drjT] y v F^(x + r)y) 



(3.2) 



(3.3) 



The parallel transport implicit in the correlator on the l.h.s. of ( |1.3| ) is taken along a V- 
shaped path consisting of two straight lines from x to xq and from xq to x', respectively. 
The path is of course not the straight line connecting x to x' which is normally used in 
lattice calculations [|l3j. One has to make the assumption that the correlation length 
is independent of the position x somewhere inside the loop jl2j.PI 



The average over the gauge field background <C .. ^> in (|3.1|) can be evaluated 
using the lowest (quadratic) term in the cumulant expansion. This is the "stochastic 
vacuum" model assumption 



< e/ A >^ exp(-l/2 f [ < AA » 



(3.4) 



The correlator in the exponential follows from eq. (|1.3| ) with (|3.3|) . After rescaling 
t = Tf, y = T x l 2 y and taking only the form factor D in ( |1.3| ) (k = 1), we obtain for 
the correlated part in (|3.1| ) 



rgauge(^ ) < 93 



xtr exp 



F 2 >) = -~ f ^T-^ J d d x J [Dy] exp{- £ df V - - m*T} 
< glF a F a > ll c 



d(d - 1) (iV c 2 - 1) Jo 



dr}dr}'(r}<r]')(-(y ■ y'f + y ■ y'y ■ y'^D 



T / dfdf' 



(rjy - rj'y'fT 



+A(l-d)l L D 



(3.5) 



+ J dr]7](y x y — y x y')D \ j + (primed *-> unprimed) > 



Expansion to order F 2 would reproduce the contribution from gauge bosons in the loop 
calculated in chapter 2. 



Some indication that the correlation length is path independent can be found in ref. 1 23 
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The ghost contribution r g h os t has the opposite sign and an additional factor of 2 
compared to ( [3.1|) . The trace is only over adjoint color indices, and only the first term 
in the curly brackets in ( |3.5D appears. 

The first term of ( |3.5| ) originates from the correlation of two Wegner- Wilson loop 
integrals exp(— § Adx). Using the nonabelian Stokes theorem one obtains area integrals 
and the correlator ( |1.3j ) between two F's for large (d 3> a) loops leads to one common 
area integral. This is a nice and well-known result of the stochastic vacuum model 



11, 12]. The first term in (|3.5|) is just this area integral. For large areas it reduces to 



a A with the (rescaled) string tension 

1 iV r°° ( z 2 \ 2N 2 

and A = A/T where A is the area of the (minimal) surface whose boundary is given 
by the closed path y(r). In eq. (|3.6[) we wrote cx a dj and o"f un d for the confining string 
tension in the adjoint and fundamental representation, respectively. 

The second term in the curly brackets in (|3.5| ) arises from the correlation of two 
explicit factors of in the exponential of eq. (|3.1|) . It is related to a gauge boson 
spin interaction which is short-ranged. Its range is of the order of the correlation length 
a. This interaction owes its existence to the nonminimal "paramagnetic" coupling of 
the gauge boson fluctuations a M (x) to the background field, i.e. to the F^-term in the 
kinetic operator K^ u . The third term in ( |3.5| ) is due the interference of the confinement 
and spin effects. 

To start, let us first disregard the "diamagnetic" interactions coming from the D 2 - 
term in and let us focus on the spin interaction. For large loops (d ^> a) the 
r, r'-correlations are well approximated by nearest neighbour interactions along the 
closed path. Instead of doing the above world-line path integral one can change to 
the ordinary QFT language and calculate the gluon self-energy £ in a correlated gauge 
background. This calculation is described in Appendix A in some detail. There we also 
include the covariant derivative induced parts of the interaction. Of course, the use of 



the correlator ( |1.3| ) again requires the specification of a reference point for the parallel 
transport operator and for the coordinate gauge. In the computation of we 
identify both of those points with x = (x + x')/2. 
For the spin-induced part of £ we obtain 

E(x, x')% = S F (z 2 = {x - x') 2 )5,J ab (3.7) 

with 

where H{^-) is a massive (m 2 = g 2 (p 2 /4) scalar propagator function. Fourier transfor- 
mation leads to 

S F (P 2 ) = +<S-F>>±fd t ,f 



pq [p — q) 1 + m z 
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where we have inserted the physical values d = 3 and N c = 2 now. The contribution 
Sf(z 2 ) is given by the last term of the expression (|A.31 ) for the complete function 
S(z 2 ). 



Eq. (|3.8| ) shows that the range of the spin interaction is determined by the fall-off 
behavior of the form factor D. Furthermore it is important to observe that Sp(z 2 ) 
is a positive function. The self-energy £ is defined in such a way that the inverse 
dressed propagator is G^ -1 — £ where is the free one (see eq. ( |A.6|) ). Hence 
Sf > implies a negative ( "tachyonic" ) contribution to the effective squared mass 
of the gauge boson. This means that the interaction with the stochastic background 
destabilizes the "empty" vacuum where the gauge boson fluctuations have a vanishing 
expectation value: They have the tendency to condense which then gives rise to a 
nonzero vacuum expectation value of F 2 . 

This conclusion does not change if we retain the interactions coming from the D 2 - 
term in K^ u , i.e. the x^A^(x)-term in the world line path integral (|3.1| ) in a coordinate 
gauge with x = (xi + x 2 )/2. The complete result is given in eq. ( |A.35| ). This delivers 
part of the gauge interaction in D 2 into £ whereas the main part constitutes the 
confining force summarized in the area law term discussed below. In eq. ([A. 38 



we 



give the complete result for the mass shift Am 2 of the gauge field fluctuations. The 
spin interaction contributes the "1" in the curly brackets of eq. ( |A.38[) , whereas the 



D 2 -interactions produce the terms proportional to 8i and 82- In view of ( |A.39|) the 
latter are seen to be negligible at our present level of accuracy. 



Inserting the self-energy (|3.9| ) in a gauge boson loop leads to an effective potential 



v ** ~ 4 / (S^ in[ °° 2 +m2 ~ ^°° 2))] (3 - io) 

Here we have multiplied by 2 because of the two different ways to get a chain of 
neighbours. This has to be corrected for the zeroth and first order terms in Sf- Sf{p 2 ) 
is positive and contributes negatively in the propagator. Of course ( |3.10| ) has to be 
UV-renormalized but the counter terms are identical to those of the case Sp = 0. Thus 
( ETlCl ) reads more correctly 

Vs * = 1 1 Wfi 2 Hp2 +m2 ~ ^°° 2 ' m2)) - Hp2 + m2) 

S F {p 2 ,m 2 ) m 2 2l 

The ghost loop- induced part has to be added - it does not contain Sf- 

Let us now focus on the other extreme and consider the approximation of a pure 
area law in the IR for the first term in (|3.5|) . Then it behaves as ~ exp(— a AT) for large 
AT. Furthermore, we assume that when AT becomes smaller this term is damped as 
~ exp(— cA n < g^F 2 > T n (a 2 ) n ~ 2 ) with an exponent n > 2. This damping is supposed 
to set in for AT < ca 2 where c = (9(1) is treated as a free parameter .f\ In the following 

5 Such a behavior with n = 2 was argued for in the context of the stochastic vacuum model in ref. 
p4[ . It would lead to a ml on{ (p 2 ) in ( 3 . 1 5| ) below with a rather weak fall-off in p 2 (like Sf(p 2 ))- The 



additional piece in S in appendix A due to ry-integrals has also such a behaviour. Thus a contribution 
with n — 2 whould be better associated with the T,(x,x') evaluated before. 
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we will use the interpolating ansatz with n = 3 

~ exp {- aA wW+TV (3 ' 12) 

If we include also an area law for the ghost loop and renormalize at m 2 = F 2 = 
we obtain 



( -vAT 3 2 \ 2 _ 

exp 5 mi — 1 + mi 



(3.13) 



where A is a complicated functional of the interagtion variable y(f): it is the area of 
the minimal surface spanned by the loop y. 

Substituting (4irT)~ 3 / 2 = J (2vF e_p2T * n ( P-13| ) as a procedure to trade T for a 



momentum integration, the T integration can be performed (numerically if c ^ 0) and 
we arrive at 

r~ = _(^|^ (4lP | B!/exp{ _|^ f £ } 



ln[p 2 + m 2 + m 2 onf (p 2 , A, m 2 )] - \np 2 - — 



(3.14) 



where m 2 onf is a complicated function defined by the requirement of producing the 
ln[p 2 + (mass) 2 ] form of the integrand. It is easy to convince oneself that this function 
always exists and is positive. It would be p 2 -independent and equal to oA if c = 
0. Thus we have an area-dependent, momentum- dependent positive m 2 onf related to 
confinement. It acts as an IR regulator in the same way as the Higgs mass but it is 
moment um- dep endent . 

The functional integral over y(f) cannot be performed in closed form. It could 
be performed numerically. Here we make the rather crude approximation that it is 
dominated by paths with areas peaked near a mean value < A >. The numerical 
value of < A > will be treated as a free parameter. For c = 2 and an average value 
< A >= 2.5 the function m 2 onf (p 2 , < A >, m 2 = 0) is plotted in fig. ||. 

The main contribution to m 2 oni (p 2 , < A >,m 2 ) comes from T<A>-values smaller 
than ~ 10 — 20g^ 4 which is not very large. With increasing values of p 2 even smaller 
T<A>-values dominate. At very large areas we would expect gluon loop holes in the 
area corresponding to the splitting of the confining string and virtual production of 
glueballs. In the case of static valence fields this would lead to a screening and a 
flattening of the linear rise of the 2-dimensional potential. As recently measured on 
the lattice this only happens in the case of very large distances d ~ 5 — 10g^" 2 much 
bigger than the inverse glueball mass ~ a/2 ~ (2 ■ 0.73(7f )" 1 The relation between 
the 2-dimensional static potential and our effective action is not clear but we would 
expect that screening is important only at very large areas TA. The screening effect 
corresponds to higher correlations and is not described by the stochastic vacuum model 
with Gaussian correlations. 
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The spin-induced £ and the confinement effects have opposite impact on the "mag- 
netic mass" . In (|3.5|) there is also the last term which arises from the interference of the 
two effects. It changes sign under the substitution y(f) =>- y(l — f) and thus the term 
would cancel if it is small. Leaving it aside one can combine the first (confinement) 



and second (spin interaction) term in the exponential in (|3.5| ) and add a corresponding 
expression without spin interaction for the ghosts. Thus one arrives at 



V 



X 



-aAT 3 



■ exp 



fc4)2 + T2 ' 9 



+ - < g\F 2 > T 



A- 

-aAT 3 



6 

-I — exp ^ 
2 P Wc4) 2 + T 2 



(3.15) 



This is an interesting structure but it seems to be impossible to treat the functional 
integral analytically. For our numerical estimates we replace the |/-integral by the 
integrand at A =< A >. We will come back to this point in the discussion in section 
4. 

If the spin-spin interaction has a shorter range than the confinement effects one 
could work with the integrand of ( |3.11|) as a rough estimate for momenta p > 1/a. At 
even shorter ranges (large p ^> 1/a) one might suspect that (|3.11| ) is not appropriate 
any longer: The correlation via the form factor D in ( |1.3|) can be neglected and we 
have the case of a constant uncorrelated background. If we assume this background to 
be pseudoabelian, it can be accounted for by an Euler-Heisenberg type formula. In the 
case of d = 3 we evaluate [21, 17 1 the gauge field and ghost contributions as^] 



V 



EH 



9s\B\ 



8tt 3 / 2 jo 
4sinh(c/ 3 |S|T 



J™ dTT- 3 ' 2 e- m2T { 



9-6 



sinh(s 3 |£|T) gz\B\T 
1 f°° dT_ a/0 m2 



+ 



2 8vr 3 / 2 



— T" 3/2 e" 
T 



m 2 T 



(3.16) 



After renormalization, the last term is just the well-known cubic term — j^(m 2 ) 3 / 2 ~ (p J 
The first integral in ( 3.16 ) is UV and IR convergent. For m 2 < g$\B\ we have an 
imaginary contribution due to a cut in (m 2 — g^\B\) 1 / 2 signalizing instability of the 
background field. This seems to be similar to the case Sf{p 2 ,tti 2 ) > m 2 in eq. ( |3.10|) . 
As is clear from the above, ( |3.16| ) is only justified below some IR cut-off T ~ a 2 /u 2 , 
i. e. for momenta p > p = u/a, where v = 0(1). One can again replace the T 
integration by a p integration using the T~ 3 / 2 substitution trick. (|3~T6| ) can be brought 
to the form ( |3 . 1 1 ) with a new polarization function Sb{p 2 ,m 2 ) instead of Sf. Sb is 
compared to Sf in fig. |j. For p 2 ~ 2/ a indeed both quantitites converge to each 
other. This is as expected since for these momenta only the lowest order in Sf ~ F 2 
contributes effectively and the nonlocality of the condensate does not play a role any 
more. 



^Renormalization in the first of refs. ||2lf after eliminating the unstable mode is not correct. 
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Figure 4: The Euler-Heisenberg Sb{p 2 ,m 2 ) compared to SV(p 2 ,m 2 ) for m 2 = (all in 
g\ units) 

4 Evaluation and Discussion 

As outlined in the previous chapter, even the 1-loop effective action generated by gauge 
bosons and ghosts propagating in the nonperturbative vacuum which is modelled by a 
correlated gauge field and a constant Higgs field background is a complicated expression 
which cannot be evaluated analytically without drastic approximations. From our 
experience with perturbative calculations of the thermal electroweak effective potential 
we expect that also (perturbative) 2-loop contributions are necessary in order to arrive 
at reliable (gauge-fixing and renormalization point insensitive) results. It also might 
turn out necessary to include higher cumulants in a description of the long-wave length 
modes by means of the stochastic vacuum model. 

If we had such an improved effective potential at our disposal, we then could actually 
determine the "local" gauge condensate < g 2 F 2 > as a function of ip by minimizing 
the effective potential in F 2 and in (p 2 oc m 2 . Without a detailed evaluation of the 
pertinent world-line path integrals (using lattice methods, for instance) we only can 
advocate a qualitative picture: Due to spin forces the state with < F 2 >= is unstable 
while the true vacuum is stabilized at some value < F 2 >^ due to confinement forces. 

In the hot electroweak phase (<p 2 ~ 0), we learned from lattice studies that the con- 
fining string tension is approximately that of pure S77(2)-YM theory. In the stochastic 
vacuum model there is the relation (|3.6| ) between < g%F 2 > and the string tension a. 
Thus we may deduce the value of the condensate from the lattice value of a rather 
than by minimizing the effective potential. 

For m 2 ^ the destabilizing effect of the 1-loop F 2 term diminishes (see eq. ( j2.20| )) 
and at some m 2 , the < F 2 >= vacuum stabilizes. Clearly < F 2 > depends on m 2 . 
Again m 2 , could in principle be read off from a balance between tree and loop F 2 terms, 
but as we said before, in praxi a 1-loop calculation is not very reliable. We thus will 
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Figure 5: m 2 on{ {p 2 ^A = 2.5, m 2 ) and Sp{j) 2 ,m 2 ) plotted a) for m? 
we used form factor 



and b) p 2 = 0; 



parametrize the m 2 -dependence of < F 2 > by 



< F 2 >(m 2 ) =< F 2 >(0)/i(m 2 



(4.1) 



where the function h satisfies h(0) = 1, h(mfy = 0, h'(mfy = 0. We choose it in 
the form 

h(m 2 ) 



,tc m 



cos 



24, 



2 mo' 



(4.2) 



"Vfund(ra 2 



0). 



The numerical values 
For the 
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4a 1 so that F 2 (m 2 ) has a mean "decay 



for m < m , and fixF] < gfF 2 > (0) 

o"fund(^ 2 = 0) ~ 0.13(?3, 1/a ~ 0.73(^1 are known from lattice studies 
new parameter m in (fl2|) we choose m 
width" 2/o. 

In general the correlation length a also will depend on (p 2 oc m 2 . For large m 2 
the perturbative part of the form factor D (with scale 1/m) which is not included in 
our analysis (because it is represented by the usual Feynman graphs) dominates. A 
dependence a(m 2 ) should come out from a self-consistency equation for the stochastic 
vacuum |33| , but this is beyond the scope of this paper. In the practical evaluation 
of our formulas a is always accompanied by a factor < g 2 F 2 > (m 2 ) which we had 
to parametrize anyway in lack of a fully dynamical treatment. Thus we keep the 
correlation length constant .0 

As outlined in section 3 and appendix A the gauge boson vacuum polarization 
E contains a destabilizing spin force effect Sf- In the case of a nearest neighbour 
interaction in a sufficiently extended loop we obtain a negative IR magnetic (mass) 2 . 

7 In all the numerical evaluations we present in section 4 and appendix A we used the form factor 
L>W of eq. ([Hi]). 



A simple substitution — 



(m 2 



a 2 (0) 



together with (ff.l|)/(|4~^) turned out to give a very 



rapid fall-off of m^ onf (p , m 2 ) and Sf(p 1 , m ) in m . 



18 



S also contains a confining part. The latter we have summarized into a (modified) area 
law exponential in the loop path integral ( [3.13j ). It leads to a positive IR magnetic 
(mass) 2 , m 2 onf (p 2 ). Such a positive (mass) 2 in the IR region acts as an IR cutoff. 
Indeed, the soft gauge quanta are already taken into account by the stochastic vacuum 
background. They should not circulate in the gauge boson loop of the effective action. 

There seems to be a paradoxon: On one hand, in order to destabilize < F 2 >= 0, 
the spin forces have to dominate the F 2 radiative correction calculated in section 2. 
On the other hand confinement forces should stabilize F 2 at some value < F 2 >^ 
and there these forces should dominate, i.e. the singularity in p 2 shifted towards 
positive p 2 by the Sp term in (|3.10|) , (|3.11|) should not be reached because now the 



positive m 2 on{ dominates. However, this is not a contradiction if Sp(p 2 ) and rn 2 orii {p 2 ) 
have different p 2 dependences and dominate at intermediate and small (p <~ -) loop 
momenta, respectively. For suitable values of < A >, this can be indeed the case as 
we see from fig. [5| Here we have chosen < A >= 2.5, c = 2 in (|3.12|) and m 2 = 0. 

It is important to note that also Sp is changed (diminished!) by the effective IR cut- 
off through m 2 onf ((p±g) 2 ) entering (|3.9[). If m 2 onf is too large, Sp is strongly suppressed 
and never dominates - this would exclude the destabilization of < F 2 >= 0! 

The size of the m 2 on{ (p 2 = 0) chosen above is of the same order as the gap mass mNjL 



of refs. |25] who give a pole mass ttt-njl = 0.28^3 for "the magnetic mass" in a 1-loop 



Nambu-Jona-Lasinio-type equation for the gluon propagator. It is also in the range of 



the Landau gauge gluon propagator mass, m L = 0.35g 2 in the lattice evaluation of [30 



and of the magnetic mass obtained from gap equations for pure YM theory [j3~2jl ,|35"l. 
Extrapolating our rn 2 oni {p 2 ) to the "pole-mass" (at negative p 2 = —ml onf (p 2 )) would 
increase its value according to fig. || 

We can bring together the negative —Sf(p 2 ) and the positive m 2 on{ (p 2 ) in one ex- 
pression in order to obtain a rough approximation to the gauge particle contribution to 
the 1-loop effective potential. We also add the ghost part which contains the m 2 onf (p 2 )- 
modification only: 

V? = lf^{2ln[p 2 + m 2 + m 2 COQ{ (p 2 , m 2 )-Sp(p 2 , m 2 )} 



2\ ™2 



i r 2 , 2 , 2/2 2m , $f(P , m m i 2l 

In [p + m + m conf (p , m ) + ^— 5— 5 — ^5 ^- tt - hip \ 

L 1 p 2 + m 2 + m 2 onf (p 2 ,m 2 ) p 2 > 

'ln[p +m +m conf (p,m)J - hip > (4.3) 



2 Jo 2tt 2 [ L " conlv " ' n p 2 

with the modified Sp of (|3~^ 



c/2 2^ <g 2 3 FF>(m 2 ) r°° 2 1 

^ (p ' m) = 9^ Jo dqq m 

in f i P + ^:^: m rii P + 9 f 2 ' m2 2 | ) D iq 2 ) (4.4) 
V VP - Q) +m 2 + m 2 onf((P ~ 1), m ) J 

As we argued before (fig. §) also the Euler-Heisenberg region p > - is properly in- 
cluded. 
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Figure 6: Vf(<f), V/ 1 \ip) compared to the perturbative tpi term 



A similar 1-loop potential due to scalars (Higgs plus Goldstone bosons) in the loop 
has the form 



V* 



+3 



dpp 2 

2 1 _2 



ln(p + m 3 + 3\tp + m conf(fund) (p , m )) - hip - 



2 m| + 3A<^ 2 



ln(p + m 3 + Ay? + m + m conf(fund) (p , m )) - hip 



m 2 , + Ay? 2 + m 2 

^2 



(4.5) 



m conf(fund) * s related to m 2 on{ via the relation between the string tension in the adjoint 



and fundamental representation ( |3.6| ) as 



2 / 2 2\ 

m conf(fund) IP 5 m J 



-m 2 con{ {p 2 ,m 2 



(4.6) 



In the evaluation of (|4.5|) we vary m 2 starting from positive values (high temperature) 
and avoid the unstable part with negative m 2 where a treatment in the framework of 
coarse grained actions seems to be indispensible. 

The 1-loop potentials (4.3) and (|4.5| ) are the main tool for our investigation of the 



phase structure of the hot electroweak theory. They are plotted in fig. 6 for A = 0.1 lg 2 
and m 2 = — 0.0041g| together with the perturbative y3 3 -term. They have to be added 
to the tree Higgs potential 

V^dl) 



(4.7) 



written with the dimensionless quantities (|L2 ) and the dimensionless field y^ = 
<p 2 /g 2 {m 2 = g 2 ip 2 /4 = pI^dl/^)- Vpp e of eq. Q2.33] ) turns out to be small numerically. 

As we can see in fig. 7, starting at small x the potential ( |4.7|) signals the usual first- 
order phase transition, i.e. two degenerate minima at the appropriate y (temperature 
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Figure 7: Fading away of the first order phase transition with increasing x = -4, where 
x x = 0.06, x 2 = 0.08 and x 3 = 0.11 



T c ) and in between a bump with a maximum. Increasing x the bump between the 
minima and the "size of the broken phase" <p m i n becomes smaller. In the case A = 2.5, 
c = 2 and tuq = 4a _1 the bump vanishes at x ~ 0.11. This is in agreement with lattice 
calculations: we have a crossover. For this result the Sf spin part is very important .0 
In the work of Buchmiiller and Philipsen [^| the crossover point was determined 



by comparing the "magnetic mass" m^jL of the gauge bosons in the hot phase with the 
(perturbative) gauge boson mass in the broken phase in the unitary gauge. Since our 
"magnetic mass" is momentum dependent this is not a clear prescription in our case. 
Trying a similar reasoning we could compare m 2 onf (p 2 = 0, m 2 = 0) — £(jo 2 = 0, m 2 = 0) 
with the perturbative gauge boson mass in the "broken" phase. 

At the endpoint of the first-order phase transition line in the x — y plane ("crossover 
point") we expect a second-order phase transition and thus conformal symmetry. The 
latter leads to the requirement that the dimensionful parameters m 2 and A e ff mul- 
tiplying the <f) 2 - and 4 - terms of the effective 3-dimensional scalar action should be 
zero at the crossover point. At the present level of accuracy we can be gentle about 
anomalous dimensions which may indeed be rather small [26fl. Taylor expanding the 



potential V of eq. ( |4.7| ) allows us to read off the coupling A e g and with the parameters 
chosen in fig. 7 we obtain 

x cro = 0.11 (4.8) 

It is an interesting question whether the second-order phase transition is in the uni- 
versality class of scalar 4 theory and the Ising model [|26|].[^] Of course our approach 
is not quite adequate to give a satisfactory answer to this. 



9 For the parameter set used in fig. [?] "^onf — Sf is of the order of rri^ ]L obtained in ref. |2q| . The 
Higgs part of the potential just contains l^conf which is a ls° of that order whereas the ghost part 



TO co„f 

10 This was reported recently as a lattice result E7 
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The table below contains the contribution of Vf and V^ 1 to the 1-loop ip 4 term 
and the effective coupling at the x values of fig. [7|. We also give the surface tension 
in the thin wall approximation cr s = ij2V((p)dip for the full and the perturbative 
potential. We get the qualitative features of the comparison of perturbative and lattice 
results for a* [El. 



X 


5x 9 


5x H 


Xefi 






0.06 


-0.063 


-0.039 


-0.041 


0.016 


0.013 


0.08 


-0.063 


-0.043 


-0.026 


0.004 


0.007 


0.11 


-0.063 


-0.051 








0.004 



Our most important calculation of the F 2 term in chapter 2 only contained the 
correlation length a which we took from lattice data; there is also some uncertainty 
in this term because of its gauge dependence which should be reduced by a 2-loop 
calculation. The construction of an effective potential in 1-loop order contained many 
more parameters: < F 2 > (0) and the function h(m 2 ) containing tUq in (4.1) could 
in principle be obtained by studying the nontrivial minimum of the effective potential 
in F 2 ] c and n = 3 in fl3.13|) are related to the modification of the area law at small 
areas. < A > could be determined in a numerical study of the path integral ( |3.15| ). 
Given our ignorance about the last three parameters c, n and < A > we just fitted a 
potential in order to reproduce the crossover at the x- value given in lattice calculations. 
We then did not consistently determine < F 2 > (0) and h(m 2 ) from the potential but 
fixed the former by lattice data on the string tension and used a plausible value for m§. 
This point requires more exhaustive numerical studies, not to speak about a 2-loop 
calculation again. 

We have described the IR effects of the gauge interactions by use of the stochastic 
vacuum correlation ( |1.3|) . This is quite different from a renormalization group treat- 
ment like in ref. || where differential renormalization group equations (at least in 
principle) lead to an effective IR action starting from the usual tree action in the UV. 
The stochastic vacuum prescription is a pure IR concept and - in a sense - complemen- 
tary to the renormalization group approach. The latter is applied most easily if the 
relevant momentum scale is not too far below the UV cutoff where the initial condi- 
tion (classical action) is specified. In this situation simple truncations of the space of 
actions are quite successful since the effective (average) action does not deviate very 
strongly from the classical one. Once the variable IR cutoff enters the deep IR-regime 
a more exhaustive parametrization of the space of actions becomes indispensable. For 
pure Yang-Mills theory, a first step in this direction has been done in ref. [p4}| , but 
a complete ab initio treatment of the IR physics is necessarily quite involved and not 
available yet. 

The stochastic vacuum model, on the other hand, makes a phenomenological ansatz 
for the low momentum modes of the gauge field. In principle this effective IR dynamics 
should follow from the renormalization group analysis but clearly the actual proof is 
a formidable task. Accepting the stochastic vacuum as a starting point one is in the 
opposite situation as in the renormalization group framework where one is given an 
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effective average action at a certain scale k. In the former case, only the modes with very 
small momenta are described by the correlators used and all other modes with larger 
momenta, have to be dealt with ("integrated out") explicitly, e. g. by perturbation 
theory in the background of the low-momentum modes. This is what we did in sec. 2, 
for instance, when we stochastically averaged the one-loop effective action. Conversely, 
in the renormalization group case the modes with momenta smaller than k are the ones 
which still have to be integrated out. 

It would also be very interesting to see if the stochastic vacuum model can be 
derived from a self-consistent set of Dyson-Schwinger equations as proposed in ref. 

n. 

Our approach is also quite different from the Nambu-Jona-Lasinio-type equations 
for the "magnetic mass" which are proposed in [25 1. In the latter approach, this mag- 
netic mass seems to have a similar meaning as the confinement mass m 2 oni introduced 
in our discussion. We argued that this is an effective (mass) 2 , which is only important 
for very small momenta and we have an additional negative (mass) 2 , —Sp, related to 
the gauge boson spin-spin interaction which is in a delicate balance with the confine- 
ment mass. This avoids the appearence of a tachyonic pole in the gluon propagator 
related to the well known IR instability of YM theory. A higher (two) loop calcula- 
tion in the nonperturbative gauge field background would be desirable (though very 
demanding) both in order to reduce the gauge dependence and in order to treat the IR 
instability properly. Indeed this also seems to be an important point in the discussion 



of gauge invariant gap equations [32, 25 1 in a recent paper of Cornwall [35|. It would 



be interesting to include a gauge background in the gap equations of ref. fl25 |. 

In our approach the nonperturbative dynamics mainly enters through the correlated 
gauge field background like in pure gauge theory but with a further Higgs field param- 
eter m 2 oc ip 2 . We have been using lattice results for the values of a and cr fund (m = 0). 
It might be possible to explore the m 2 dependence of these quantities on the lattice. 

In finite temperature QCD (and in ordinary QCD) we have no Higgs expectation 
value Lf 2 at our disposal which could act as an IR regulator and the 4-dimensional gauge 



coupling is big. Thus the status of naive perturbation theory is even much worse |28 



Still an IR-improved perturbative picture using a correlated gauge field background 
might be possible |29j and the methods introduced in the present paper may be still 
relevant. This deserves further investigation. 
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Appendix A 

In this appendix we calculate the stochastic average of the gauge boson propagator 
G = K^ 1 with the kinetic operator K a ^ v given in eq. ( |2.16|) . We shall determine G to 
first order in the condensate < g 2 F 2 >. 

To start with, we decompose K pv according to 

= K$ + K$ + K$ (A.l) 

where the superscripts denote the powers of F^ y which are contained in the respective 
terms once the gauge field 

A n( x ) = f dr] r)x v F VfM {r)x) (A.2) 

J 

is inserted. Thus we employ the Fock-Schwinger gauge centered at x p = 

rf£ = (f + m 2 ) V 

K v$ = -g(PaA a + A a p a )5 pu + 2igF (iV 

K$ = g 2 A a AJ llv (A.3) 

with p p = — id [j, and with A^ and F pu considered as matrices in the adjoint represen- 
tation of the gauge group. 

To order F 2 , the propagator then reads 

G = -G®kMg® - G®kV>G® 

+GWkWgWkWG<® + 0(F 3 ) (A.4) 

where G (0) = [K^]' 1 is the free propagator. 

Using (|1.3p together with <C 1 ^>= 1 and <C F^ ^>= the stochastic average 
G G ^> is seen to be 

G = G (0) + G (0) SG (0) + 0(F 3 ) (A.5) 
G' 1 = G ( - 0) ~ 1 -Z + 0(F 3 ), (A.6) 



i.e., 



with the "mass operator" 

E =< K^G^kU » - < > (A.7) 

The quadratic averages in ( |A.7| ) have to be performed using the correlator (|1.3| ). For 
simplicity we restrict ourselves to the case k = 1. Furthermore, we choose the reference 
point Xo = 0, which coincides with the point at which the Fock-Schwinger gauge is 
centered. As a consequence, the path-ordered exponential in eq. ( |1.4| ) equals unity 
and it follows that 

« g 2 F« u (x)F b pa (y) »= ^ ab 

x [5^5 ua - 5^5 up ]D((x - yf/a 2 ) (A.8) 
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with 

<g 2 F 2 > ,„ , 

w s i(i-i)(jv?-i) (A ' 9) 

It is convenient to write 

S = S! + S 2 (A.10) 

Si =< Q >, Q = K^G^KV (A.ll) 

S 2 = - < fT (2) > (A.12) 

and to further decompose 

Q = Qaa + <5ff + Qaf (A. 13) 



with 



[Q^W = g 2 (p a A a + A a p a )G^(p)(p^ + A^) (A.14) 
[QffU = (2w?^«G%(j))Fp u (A.15) 



[QafU = -2ig 2 {{p a A a + A a p a )G { ^(p)Fp u 

+F lia G<$(p)(p f) Af ) + Appf,)} (A.16) 

Working in a position-space representation, quantities such as Ex or Q are kernels with 
respect to the space-time coordinates and matrices with respect to Lorentz and adjoint 
group indices. Before we can evaluate the average 

Ei(x 2 ,xi)J =« Q(x 2 ,x 1 )^ » (A.17) 

we have to determine the position-space matrix elements 

QKzi)£ =< x 2 \Q^{x,p)\ Xl > (A.18) 

for a fixed background field. For the contribution coming from Qaa, say, one finds 

Q A A(x 2 ,x 1 )f u = -g 2 (T c T d ) a % 1/ 

x '4A c p (x 2 )Ai(x l )d p d a G < -°\x 2 -x l ) 
- 2A c p (x 2 )dA d (x 1 )d p G {0 \x 2 - Xl ) 
+ 2dA c (x 2 )Aj(x 1 )d a G ( -°\x 2 - Xl ) 
-dA c (x 2 ) dA d (xi ) G (0) {x 2 - x x )] (A.19) 

with dA a = d^A^ and < x 2 |G , | i °2l a; i > = $npG^\x 2 — x\). The next step is to express 
all A^s in terms of F^s by using ( |A.2j ), and to apply ( |A.8|) to the resulting <C FF 3> 
correlators. In particular one needs 



< g'Al( X2 )A*( Xl ) >= u)tir[(xi ■ x 2 )6 ia , - x lp x 2u ] 

x f dr} x \ dr} 2 r} x r} 2 D{b 2 1 a 2 ) 
Jo Jo 

« g 2 A c p (x 2 )dA d (x 1 ) »= -2^5 cd a^ 2 [(x 1 • x 2 )x 2 , - (x 2 ) 2 x lfl ] 
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x / d Vl / dr l2 {r ll r l2 ) 2 D'{5 2 /a 2 ) 
Jo Jo 

< 5 2 &4 c (x 2 )&4 d (:ri) >= -2ujb cd a~ 2 

x\(d-l)( Xl -x 2 ) f d Vl j 1 dr l2 {r ll r l2 ) 2 D'{5 2 /a 2 ) 



+2a' 2 v l2 £ d Vl £ dr ]2 {r h r ]2 YD"(8 2 /a 2 )} (A.20) 

with the abbreviations 

<5 M = V2(x 2 )^ - ?7i(xi) M (A. 21) 

and 

v 12 = (x 1 ) 2 (x 2 ) 2 -(x 1 -x 2 ) 2 (A.22) 

The derivatives acting on in eq. ( |A.19| ) can be easily removed by noting that the 
free propagator has the structure 

G®(x 2 -x x ) = H Q(x 2 - xx) 2 ) (A.23) 
for some function H. In the massless case, for instance, 

H(v) = ^(2vr)- d / 2 r(d/2 - 1) v~ d / 2+ \ d>2 (A.24) 
After some calculation this leads to 

« Q A A(x 2 ,x 1 )f I/ »= -2uN c 5 a \ u 



x |2[(d - ■ x 2 )H' + v 12 H") £ d^di] 2 i] x i] 2 D{b 2 la 2 ) 

+4a- 2 v l2 H' ^ dr ]l dr ]2 {r ]l r ]2 ) 2 D'{5 2 /a 2 ) 
Jo 

+{d - l)a- 2 (x x ■ x 2 )H j dr ]l dr l2 {r] l r] 2 ) 2 D , {5 2 /a 2 ) 

Jo 

+2a~ A v 12 H jT 1 dr jl dr ]2 {r] 1 r] 2 ) 3 D"{8 2 /a 2 )\ (A.25) 

Here a prime denotes the derivative with respect to the argument, and the function H 
and its derivatives are always evaluated at the point (x 2 — Xi) 2 /2. Likewise one finds 
for the other contributions to Ex: 

< Q FF {x 2 , Xl )% > = 4(d - 1 ViV c <5 a V#((x 2 - x 1 ) 2 /2) J D((x 2 - xrf/a 2 ) 

(A.26) 

<g AF (x 2 ,x 1 )J > = ^g^uN^x^x^ - x Xv x 2v ) 



■{h'((x 2 -xrf/2) j drjrjDdxt-rjx^/a 2 ) 

+2a- 2 H((x 2 - xi) 2 /2) drirfD'Hxx - i]x 2 fja 2 ) 

Jo 

+ (xi^x 2 )} (A.27) 
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The contribution £2 is a pure contact term which vanishes for x\ 7^ x 2 . 

M^x^t = -(d-l)uN c 5 ab 5^x 2 1 5(x 2 -x 1 ) 



f 1 drndrh-niThDifa - r) 2 ) 2 x\la 2 ) (A.28) 
Jo 



The complete "mass operator" E is given by the sum of the terms in eqs. (|A.25 ), 



( |A.26| ), (|A.27|) and ( |A.28|) . We observe that for the gauge chosen, E is not translational 
invariant; it depends on x\ and x 2 separately and not only on the difference X\ — x 2 . 
Generally speaking E is a function of the two translational invariant combinations of 
X\,x 2 and the point xo at which the Fock-Schwinger gauge is based (cf. eq. ( |3.3|) ) and 
which was chosen to act also as the reference point in the <C FF ^-correlator, see eq. 
( p..3| ). Without loss of generality, one of these three points can be fixed at will. We 
have exploited this freedom in order to set Xq = 0. The general case is recovered by 
substituting x\, 2 x\ y2 — xq. We see that E is proportional to the unit matrix in color 
space. It is not, however, a unit matrix with respect to the Lorentz indices since Qaf 
has a nontrivial tensor structure involving x\ and x 2 . 

Next we consider E with a different choice of the point xq which, again, serves both 
as the reference point in the correlator and as the base point of the Fock-Schwinger 
gauge. We define xq to be in the middle of the straight line connecting x\ to x 2 . 

xo = ~{xi +x 2 ) (A.29) 

For this choice of xo, E depends on x\ and x 2 only via [x\ — x 2 ) 2 . One finds 

S(x 2 , = 5^5 ab S(( Xl - x 2 ) 2 ) (A.30) 

The function S can be determined from our first calculation by setting 

1 1 

Xl = -z, x 2 = --z 

for some vector z. Then xo = \{x\ + x 2 ) = 0, and the reference points agree. In this 
manner one is led to the comparatively simple result (z = x\ — x 2 ) 



■{z 2 H' ^-j drjidrj^^D 




dr]idr] 2 (r]ir] 2 ) 2 D' I (771 + rj 2 f 




The Qaf contribution vanishes. 

Since the operator E calculated according to the second prescription is translational 
invariant, it is meaningful to Fourier-transform it with respect to z, and to interpret 
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,m 2 ), 



Figure 8: Defining S(p 2 ,m 2 ) =< g 2 F 2 > a 2 ■ s(p 2 ,m 2 ) the full line is s(p 2 = 
the dashed and the dot-dashed lines are the contributions of the spin and the minimal 
interactions to s(p 2 = 0,m 2 ), respectively. We used the form factor 



E(p 2 = 0) as a kind of mass shift originating from the interaction of the gauge bosons 
with the stochastic background field. If we define 



s{v 2 ) 



d d ze- ipz S(z 2 ) 



(A.32) 



and use the Fourier transforms 



D P. 



d d q 
d d q 



e iqz (q 2 + m 2 )~ 1 



e iqz D{q 2 ) 



(A.33) 



it is straightforward to arrive at 

Nr. 



S(P 2 



< g 2 F 2 > 



d(N 2 - 1; 



d d q f D{q 2 ) 



-T I dri x di] 2 



V1V2 



(2vr) fl! [ (p - q) 2 + m 2 
d-2 2m 2 



(V1V2) 2 
(Vi + V2) 2+d 



(p — q) 2 + m 2 [(p — q) 2 + m 2 ] 2 
4q 2 \ 



D 



4q 2 



(Vi + V2)' 



dD 



8q 2 



:D' 



4q 2 



(rjx + V2) 2 , 
1 



(V1 + V2) 2 \(vi + V2) 2 



(p ~ q) 2 + m2 



(A.34) 



For the remaining analysis we restrict ourselves to the physically most interesting case 
in this paper, d = 3 and iV c = 2. The angular integration in (|A.34| ) can be evaluated 
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analytically. One obtains 
2 



9tt 2 



< qIf 2 > 



dqq 2 



—D(q 2 ) In + if + m2 
2pq (p — q) 2 + m 2 



fa + w 



1 , (p + q) 2 + m 2 
— In ±L — ±1 1 

pq (p — q) 2 + m 2 (p 2 + q 2 



5WI 



■D 



4g 2 



v(m + ^) 2 

, , (V1V2) 2 1 
^1^27 — ■ — tt— 
fa + ^fpq 



3D 



4q 2 



+ 



8g 2 



D' 



2\2 



Ap 2 q 2 



4q 2 



(V1+V2) 2 / (V1 + V2) 2 \(vi+V2) 2 



•In 



(p + q) 2 + rn? 
(p — q) 2 + m 2 



(A.35) 



The integral representation ( |A.35| ) is the main result of this appendix. The remaining 
g-integration has to be performed numerically. Results are given in fig. As we 

can see the first term due to the spin interaction dominates. The other terms related 
to minimal interaction have opposite sign and go into the same direction as the area 
law-induced mass square discussed in section 3. 

Indeed, instead of the choice ([A.29Q for the base point xq one could choose xq to 



be the center of mass introduced in eqs. (|3.2|) . The 77-integrals together with the r- 
integrals in the first part of this appendix then belong to area integrals filling the loop 
in the sense of the nonabelian Stokes theorem. This is similar in spirit to ref. [33j. One 
could imagine a split of these area integrals into segments related to the ^-contributions 



to ( ]A.35|) and a main area integral approximated as in section 3 by a (modified) area 
law. 

We mentioned already thatQ Am 2 = —S(p 2 = 0) should be interpreted as the mass 
shift due to the interaction with a stochastic background. This quantity turns out to 
be IR- finite even in the case when the gauge bosons are massless a priori, i.e., setting 
m = on the RHS of eq. (|A.35|) yields a finite function S(p 2 ) and in particular a finite 
value for 5(0). It is even possible to calculate this value analytically. Returning to the 



representation ( A.31|) , one obtains for arbitrary m 2 , d and A^ c : 



-Am 2 = / d d zS(z 2 ) 
•1 



71 



d/2 N 



< g 2 F 2 > 



dT(d/2)(N 2 



r 



» w \ 



dvidmvm J o dw w d/2 l wH' (-) D f + mf-^ , 
P dvidri2(viV2) 2 / o °° dw w d / 2 -\w /2a 2 )H (^j D> (fa + rj 2 ) 



+ 



w 



Aa 2 ) 



(A.36) 



11 For the form factor D^ 2 \ s(0,p 2 ) differs by roughly a factor of 2 but the relation between < g 2 F 2 > 
and crfund also changes and the result for S is close to that for 

12 The minus sign is due to the fact that the inverse dressed propagator is G^ -1 — S. 
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Assuming now the propagator function H to be of the massless form (|A.24 ), the inte- 



grals in ( |A.36|) can be simplified considerably. In this case one may use the identities 

(d>2) 



w 



d/2~l H f^\ = i 7T -d/2 TW2 _ 1) (A37) 



in order to bring ( |A.36| ) to the following explicit form: 

2/V 

A 2 z,iv c 2 . 2 T->2 ^ 

" Am = d(d-2)(iV?-l) fl <^ > 



/■oo 

•{1 - (d-2)0i -29 2 } / cM>(w) (A.38) 

Jo 



with the constants 



WoM'%!r^ = ^= 0X1417 --- (A ' 39) 



The remaining integral in (|A.38| ) slightly depends on the precise shape of the model 



function D we have chosen. For the exponential D(w) = e~ w } corresponding to eq. 
( [2.12j ) the integral J Q °° dwD(w) equals unity, for instance. In any case it is a positive 
number of order unity. 

Since in the physically interesting situation d = 3 the expression inside the curly 
brackets of ( |A.38|) is positive, we are led to the important conclusion that Am 2 is 
negative (tachyonic). This shows that the interaction with the stochastic background 
destabilizes the perturbative gauge field vacuum and drives the system towards the 
formation of a condensate. 



Appendix B 

In this appendix we explain the basic mechanism of why stochastic background fields 
lead to the formation of a condensate within a simple scalar toy model. 

Let be a real scalar Gaussian random variable, i.e. the average of an arbitrary 

functional F[J] is given by 



« F[J] »= J VJexp(-^J JQJ)F[J] (B.l) 



where 



J jnj = J d d xd d yJ{x) < x\Q\y > J{y) (B.2) 
with some positive definite operator Q. Eq. ( |B.1| ) implies that for any function <f>(x) 

« exp (± J d d xJ(x)<p(x)^j »= exp (+- J <j>D$\ (B.3) 
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where the operator D = VL~ l is positive definite, too. We assume that D (and fl) is 
translational invariant with real matrix elements < x\D\y >= D(x — y). Then the only 
nonzero connected correlation function reads 

« J(x)J(y) »= D(x-y) (B.4) 



This equation should be compared for the correlator (|1.3|) of the stochastic vacuum 
model in QCD. 

Let us consider the partition function Z[J] of a scalar field theory governed by 
the action S[4>] with an additional linear coupling of the field <fi(x) to a stochastic 
background J(x): 

Z[J] = J £>0exp l-S[(f>] - J d d xJ(x)(/)(x)\ (B.5) 

Its average is given by 

«Z[J}» = J D<j>e- S ^ « e~^ H » 

= |p0exp{-S[0] + ~|<^j (B.6) 

We see that a linear coupling of to a Gaussian random field has the effect of changing 
the classical action according to 

Sty] _> S ty] - I J d d xd d y<t>{x)D{x - y)<S>{ V ) (B.7) 



It is important to observe that the new nonlocal term in (|B.7|) is negative for any 4>(x). 
This means in particular that it supplies a negative contribution to the mass square of 
0. If S is of the Z 2 -symmetric form Sty] = J d d x ||(<9 M 0) 2 + |m 2 2 + ^A0 4 j, say, then 
the stochastic background induces a mass shift 

m 2 -> m 2 + Am 2 (B.8) 

with 

Am 2 = - J d d xD(x) = -D(p = 0) < (B.9) 

where D is the Fourier transform of D. If m 2 > in absence of the J<ft coupling, the 
theory has a ^-symmetric vacuum with the minimum of the potential at = 0. If, 
however, the coupling to the stochastic background is switched on and if |Am 2 | > m 2 , 
then m 2 + Am 2 becomes negative. Hence the modified classical potential develops two 
minima at ^ and the Z 2 -symmetry is spontaneously broken. 

In Yang-Mills theory a similar but technically less transparent mechanism is at work. 
The Lagrangian for the gauge field fluctuations a^(x) reads ^a^K^ u [A]a u + 0(af l ), and 
clearly and are analogous to and J, respectively. The coupling of a M to the 
background is more complicated than the linear J0-term, however. In particular, 
the analogue of the correlator ( |B.4|) is formulated in terms of -F^fA] rather than A^ 
itself. This problem can be overcome by employing the Fock-Schwinger gauge which 
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provides a simple formula for in terms of F^, see eq. ( |2.5| ). Both the computation 
of the effective action in the main body of the paper and the calculation of the mass 
operator in Appendix A show that at least for d = 3 the Am 2 term for a M which is 
induced by the stochastic background is negative, see eq. ( |A.38|) , for instance. For 
large values of d we find Am 2 > 0, however. 
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